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Abstract. The adiabatic A-soliton train interactions for the scalar nonlinear Schro- 

dinger (NLS) equation and its perturbed versions are well studied. Here we briefly outline 

how they can be generalized for the higher NLS-type equations and for the multicompo- 

nent NLS equations. It is shown that in all these cases the complex Toda chain plays 

Q ■ fundamental role. 

J£ • 1. Introduction. 

C/3 ■ 

In a number of applications to fiber optics communications (see e.g. [p], 0, [3| and the 

numerous references therein) it is important to analyze the iV-soliton train interaction of 

the (perturbed) nonlinear Schrodinger equation ((p)NLS) and some of its multicomponent 

versions. An important class of these equations are known to be integrable by applying 

the inverse scattering method to the generalized Zakharov-Shabat system 11: 

ty=(i^ + q{z,t)-\j)i/>{z,t,\) = 0, J=(l\\ q(x,t) = (l^Y {1) 

\^q \ where q(x,t) and J are matrices with compatible block structure. The case when u is a 

P-component vector is relevant to the vector NLS known also as the Manakov model |4j] ; 
Q\ | more general multicomponent NLS equations are obtained if u(x,t) is a generic rectan- 
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gular matrix. 



The class of higher multicomponent pNLS equations may be written as ||: 

£ : iJ Qt + 2f(A)q(x,t) = iJK[u], K[u] = (^ R ^ J , (2) 

S^ where /(A) is the dispersion law (for the NLS eq. /(A) = A 2 ), and A is the recursion 

operator 



A = A+ + A_, A ± x = l 



dX 
dx 



iq(x) I dy \ti ([J, q},X(y)}). (3) 

'±oo A 



In what follows we consider only dispersion laws /(A) = X 2 g(X) with polynomial <?(A). 

2. A-soliton trains of the higher scalar NLS equations. 

Let us start with a brief analysis of the A-soliton train interactions for the higher 
scalar (i.e., P = 1) NLS equations; by A-soliton train we mean a solution of the equation 
satisfying the initial condition 
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0) = ^4 s (x,0) (4) 



where u\ s (x,t) is the one-soliton solution of the higher NLS eq. Below following the 
Karpman-Solov'ev parametrization |j| we write: 

u k {x, t) = — — z k (x, t) = 2v k {x - £ fe (*)), i k {t) = —^-t + ffc.o, 

cosh{z k {x,t)) v k 

, / ,x 2/i fc . 2(fi k fl k — Vkfa,k) , . r / K \ 

where £&(£) characterizes the center of mass position of the fc-th soliton and u k , fi k and 
5 k characterize its amplitude, velocity and phase. Here f(X k ) = fo,k ± ifi,k is the value 
of the dispersion law for Xf = fx k ± iv k . As it is well known, L ([!]) remains isospectral if 
q(x,t) satisfies (0) with 1Z = 0. 

The adiabatic approximation means that the solitons are well separated and have 
nearly equal amplitudes and velocities. If we denote by e the overlap of two neighboring 
solitons and assume that at t = we have £fc(0) = ^ j0 and ^+i,o — £fc,o — r o then 

|A+-A+| 2 ~e<l, e ~ 16z/ 2 r e" 2 " oro , u r > 1, \v - i/ fc)0 |r « 1. (6) 

where A J = /i + zV = Y!k=i(^k + w k )/N and A^ = (A^)*. 

In the adiabatic approximation all iV solitons keep their identity; the main part of 
the energy of the train is related to the 2N discrete eigenvalues of L which coincide with 
A fc only in the limit r — > oo. Since the solitons are well separated we can describe the 
slow evolution of their parameters by deriving a dynamical system for them. Thus in 
[[J we derived the generalized Karpman-Solov'ev system (GKS) describing soliton trains 
with N > 2. After some additional approximations using (^) the GKS simplifies to the 
complex Toda chain (CTC) with N sites: 

iQh. = C ( e <J*+i-3* - e <fc-0*-i) + n[Q], (7) 

where C = Ag h , f = /(A+), g = g(X+), h = /(A+) - f(X ) - (A+ - X )f'(X ) and 

Qk = -2Mk + i{2iiotk-6 k ) + Q{t)+ikK, -± = 2if + ^ h . (8) 

at vo 

H[Q] is determined by the perturbative terms TZ[u] in the pNLS equation, see ||. For 
/(A) = A 2 we have g = 1, C = 16z/ 2 ; besides the linear function Q(t) can be replaced 
by t/NJ2 k =i $k and thus we reproduce the results for the A^-soliton trains of the scalar 
pNLS, see g B 1- 

Note that Q(t) can be adjusted as convenient. Indeed, it is a linear function of t so it 
does nor show up in the l.h.side of (|7|). At the same time it does not depend on k so it 
does not show up in the exponentials in the r.h.side of (|7|) as well. 

The importance of the CTC is based on the following facts: 

i) for lZ[u] — we have 1Z[Q] =0 and we deal with a completely integrable system. 
As a result, given the initial values of the soliton parameters, we are able to predict the 
asymptotic regime of the corresponding iV-soliton train, || ; 

ii) the CTC, in contrast with the well known real Toda chain, possesses much richer 
sets of asymptotical regimes. These include in addition bound states, mixed states, 



singular regimes etc. In || we singled out a special subclass of iV-soliton bound states in 
which the solitons move quasi-equidistantly. 

iii) the CTC is an universal model. Indeed, it describes the iV-soliton train for all 
the higher NLS equations. The only dependence on the dispersion law in ([!]) is in the 
constant C which can easily be taken away by changing qj, — > q^ + k In C. 

iv) If the coefficients in R[u] are of the order of e then the corresponding 1Z[Q] becomes 
a constant of the order of e. As a result the perturbation drives only the center of mass 
and the total phase of the iV-soliton train and does not influence the relative coordinates 
and phases, see 0. Our final remark here is that there are indications, showing that the 
CTC may describe iV-soliton trains also for the NLEE related to more complicated Lax 
operators, e.g. to ones depending quadratically on A, ||. 

3. iV-soliton trains of MNLS equations. 

Our main idea here is to show that some of the results in Section 2 can be generalized 
also to a special type of MNLS equations of the form: 

du p l<Pu p * A , |2 , Dr , 

l ltt + 2lhF + ^ ps ^ Up= p ^' (' 

s=l 

where A ps is a P x P symmetrical matrix. Some of these multicomponent equations with 
R[u] = are integrable. They can be written in the form ([|) and treated by the inverse 
scattering problem method applied to ([!]). Among them is the Manakov model [|| with 
A ps = (— l) Ps , where p s takes values and 1. In these cases the soliton solutions and 
their interactions in the generic case are well known ||, ||, [J. 

In what follows we consider somewhat more general case with A sp satisfying the 
condition Yl s =i ^sp = a 2 = const for all values of p. They are Hamiltonian models with 

^4£^(ffi'I)-|/->-i 2 w 2 )- < io > 

A particular case of such MNLS with P = 2 and A = ( „ J with (3 ^ 1 plays important 

role in nonlinear optics: it describes propagation of light in birefringent fibers, see e.g. 
[0, |^|. Although for such generic A the MNLS (|9]) is not integrable it allows a special 
type of solutions whose behavior is very close to the ones of the scalar NLS equations. 

Indeed, let us consider an iV-soliton train of (0) defined as the solution to (|^) satisfying 
the initial conditions 

N 

u(M) = ^Ai s (x,o), (li) 

fc=i 

where ul s (x,t) is the one-soliton solution of the scalar (P = 1) NLS equation and n^ is 
the polarization vector of the fc-th solitons in the train. If we put (n^ k ') s = e l,fik3 a~ 1 then 
obviously ^2 S=1 A sp \n s | 2 = 1 and the system (|S|) reduces to the scalar NLS equation. 

Now it remains to take into account that in the scalar case the iV-soliton interaction 
is described by the CTC. The important difference as to the scalar case is that in the def- 
inition of the corresponding q^ besides the phases (fk there enter also the scalar products 



of the corresponding polarization vectors [[UJ] : 

Q k+1 -Q k = 2i\+{£ k+1 -£ k )+iL- 4+i + 5 k + In fabl^A} . (12) 

As it is the CTC for the center of mass coordinates q k is not sufficient to determine the 
dynamics for the MNLS iV-soliton train. The complete treatment should include also an 
additional system of equations describing the dynamics of the polarization vectors n k . 
Work on this is in progress. 

4. Conclusions. 

The dynamics of the (perturbed) scalar higher NLS equations is described by the 
(perturbed) CTC. The only thing that depends on the dispersion law is the coefficient 
C in (0). If the perturbations coefficients in R[u] are of the order of e then 1Z[Q] is a 
constant of the order of e whose effect is to drive only the center of mass motion of the 
./V-soliton train and the total phase. The relative positions £, k+ i(t) — £*.(£) and the phase 
differences S k+ i(t) — 5 k (t) are not influenced by these perturbations. 

In the multicomponent case we showed that some special types of iV-soliton trains of 
the MNLS equations reduce to an extension of the CTC by additional system of dynamic 
equations responsible for the dynamics of the polarization vectors n k . 
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